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Abstract 

Relations between the length of a sunspot cycle and the average temperature 
in the same and the next cycle are calculated for a number of meteorological 
stations in Norway and in the North Atlantic region. No significant trend 
is found between the length of a cycle and the average temperature in the 
same cycle, but a significant negative trend is found between the length of a 
cycle and the temperature in the next cycle. This provides a tool to predict 
an average temperature decrease of at least 1.0 °C from solar cycle 23 to 24 
for the stations and areas analyzed. We find for the Norwegian local stations 
investigated that 25-56% of the temperature increase the last 150 years may 
be attributed to the Sun. For 3 North Atlantic stations we get 63-72% 
solar contribution. This points to the Atlantic currents as reinforcing a solar 
signal. 
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1. Introduction 

The question of a possible relation between solar activity and the Earth's 
climate has received considerable attention during the last 200 years. Periods 
with many sunspots and faculae correspond with periods with higher irra- 
diance in the visual spectrum and even stronger response in the ultraviolet, 
which acts on the ozone level. It is also proposed that galactic cosmic rays 
can act as cloud condensation nuclei, which may link variations in the cloud 
coverage to solar activity, since more cosmic rays penetrate the Earth's mag- 
netic field when the solar activity is low. A review of possible connections 
between the Sun and the Earth's climate is given by Gray et al. (2010). 

Based on strong correlation between the production rate of the cosmo- 
genic nucleids ^^C and ^°Be and proxies for sea ice drift. Bond et al. (2001) 
concluded that extremely weak pertubations in the Sun's energy output on 
decadal to millennial time scales generate a strong climate response in the 
North Atlantic deep water (NADW). This affects the global thermohaline 
circulation and the global climate. The possible sun-ocean-climate connec- 
tion may be detectable in temperature series from the North Atlantic re- 
gion. Since the ocean with its large heat capacity can store and transport 
huge amounts of heat, a time lag between solar activity and air tempera- 
ture increase is expected. An observed time lag gives us an opportunity for 
forecasting, which is the rationale for the present investigation. 

Comparing sunspot numbers with the Northern Hemisphere land temper- 
ature anomaly, Priis-Christensen & Lassen (1991) noticed a similar behavior 
of temperature and sunspot numbers from 1861 to 1990, but it seemed that 
the sunspot number R appeared to lag the temperature anomaly. They found 
a much better correlation between the solar cycle length (SCL) and the tem- 
perature anomaly. In their study they used a smoothed mean value for the 
SCL with 5 solar cycles weighted 1-2-2-2-1. They correlated the temperature 
during the central sunspot cycle of the filter with this smoothed weighted 
mean value for SCL. The reason for choosing this type of filter was that it 
has traditionally been used to describe long time trends in solar activity. 
However, it is surprising that the temperature was not smoothed the same 
way. In a follow up paper Reichel et al. (2001) concluded that the right cause- 
and-effect ordering, in the sense of Granger causality, is present between the 
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Figure 1: Length of solar cyeles (inverted) 1700-2009. The last point refers to SC23 which 
is 12.2 years long. The gradual decrease in solar cycle length 1850-2000 is indicated with 
a straight line. 

smoothed SCL and the cycle mean temperature anomaly for the Northern 
Hemisphere land air temperature in the twentieth century at the 99% signifi- 
cance level. This suggests that there may exist a physical mechanism linking 
solar activity to climate variations. 

The length of a solar cycle is determined as the time from the appearance 
of the first spot in a cycle at high solar latitude, to the disappearance of the 
last spot in the same cycle near the solar equator. However, before the last 
spot in a cycle disappears, the first spot in the next cycle appears at high 
latitude, and there is normally a two years overlap. The time of the minimum 
is defined as the central time of overlap between the two cycles (Waldmeier, 
1939), and the length of a cycle can be measured between successive minima 
or maxima. 

A recent description of how the time of minimum is calculated is given 
by NGDC (2011): "When observations permit, a date selected as either a 
cycle minimum or maximum is based in part on an average of the times 
extremes are reached in the monthly mean sunspot number, in the smoothed 
monthly mean sunspot number, and in the monthly mean number of spot 
groups alone. Two more measures arc used at time of sunspot minimum: 
the number of spotless days and the frequency of occurrence of old and new 
cycle spot groups." 

It was for a long time thought that the appearance of a solar cycle was 
a random event, which means for each cycle length and amplitude were in- 
dependent of the previous. However, Dicke (1978) showed that an internal 
chronometer has to exist inside the Sun, which after a number of short cycles, 
reset the cycle length so the average length of 11.2 years is kept. Richards et 
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al. (2009) analyzed the length of cycles 1610-2000 using median trace analy- 
ses of the cycle lengths and power spectrum analyses of the O — C residuals 
of the dates of sunspot maxima and minima. They identified a period of 
188 ± 38 yr. They also found a correspondence between long cycles and 
minima of number of spots. Their study suggests that the length of sunspot 
cycles should increase gradually over the next ~ 75 yr. accompanied by a 
gradual decrease in the number of sunspots. 

An autocorrelation study by Solanki et al. (2002) showed that the length 
of a solar cycle is a good predictor for the maximum sunspot number in 
the next cycle, in the sense that short cycles predict high Rmax and long 
cycles predict small Rmax- They explain this with the solar dynamo having 
a memory of the previous cycle's length. 

Assuming a relation between the sunspot number and global tempera- 
ture, the secular periodic change of SCL may then correlate with the global 
temperature, and as long as we are on the ascending (or descending) branches 
of the 188 yr period, we may predict a warmer (or cooler) climate. 

It was also demonstrated (Friis-Christensen & Lassen, 1992; Hoyt and 
Schatten, 1993; Lassen & Friis-Christensen, 1995) that the correlation be- 
tween SCL and climate probably has been in operation for centuries. A 
statistical study of 69 tree rings sets, covering more than 594 years, and SCL 
demonstrated that wider tree-rings (better growth conditions) were associ- 
ated with shorter sunspot cycles (Zhou and Butler, 1998). 

The relation between the smoothed SCL and temperature worked well as 
long as SCL decreased as shown in figure 1. But when the short cycle SC22 
was finished ThejU & Lassen (2000) reported a developing inconsistency. In 
order to explain the high temperatures at the turn of the millennium, the not 
yet finished SC23 had to be shorter than 8 years, which was very unlikely, 
since there had never been observed two such short cycles in a row (see fig. 
1). They concluded that the type of solar forcing described with this SCL 
model had ceased to dominate the temperature change. Since the final length 
of SC23 became 12.2 yrs, the discrepancy became even bigger. 

A more rigorous analysis was done by ThejU (2009), this time with a filter 
1-2-1 also used for temperature. In particular he investigated if the residuals 
from the relation between Northern Hemisphere (NH) land (ti2i) and SCL121 
were independent and identically distributed. He also investigated possible 
lags between SCL121 and ^121, and found that time lags of 6 or 12 years were 
necessary to obtain zero serial correlation in residuals in the data set used 
by Friis-Christensen & Lassen (1991). For updated NH-land series from the 
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HadCRU Centre, serial correlation in residuals could not be removed (ThejU, 
2009). 

In addition to the relation between solar cycle length and the amplitude 
of the next Rmax, it is reasonable to expect a time lag for the locations 
investigated, since heat from the Sun, amplified by various mechanisms, is 
stored in the ocean mainly near the Equator, and transported into the North 
Atlantic by the Gulf Stream to the coasts of Northern Europe. An example 
of time lags along the Norwegian coast is an advective delay between the 
Faroe-Shetland Channel and the Barents Sea of about 2 years determined 
from sea temperature measurements (Yndestad et al., 2008). 

Formation of NADW represents transfer of upper level water to large 
depths. The water is transported and spread throughout the Atlantic and 
exported to the Indian and Pacific oceans before updwelling in Antarcic wa- 
ters. The return flow of warm water from the Pacific through the Indian 
ocean and the Caribbean to the North Atlantic, a distance of 40000 km, 
takes from 13 to 130 years (Gordon, 1986). There appears to be solar "fin- 
gerprints" that can be detected in climate time series in other parts of the 
world with each series having a unique time lag between the solar signal and 
the hydro-climatic response. Perry (2007) reports that a solar signal com- 
posed of geomagnetic aa-index and total solar irradiance (TSI) is detected 
with various lags from zero years (Indian Ocean) to 34 years (Mississippi river 
flow) and 70 years (Labrador Sea ice). Mehl et al. (2009) have shown that two 
mechanisms: the top-down stratospheric response of ozone to fluctuations of 
shortwave solar forcing and the bottom-up coupled ocean-atmospheric sur- 
face response, acting together, can amplify a solar cyclical pulse with a factor 
4 or more. Since our stations are located near or in the North Atlantic, solar 
signals in climatic time series may arrive with delays of decades. If we can 
detect a solar signal and measure the delay for individual regions, we may 
have a method for future climate predictions. 

Recognizing that averaged temperature series from different meteorologi- 
cal stations of variable quality and changing locations may contain errors, and 
partially unknown phenomena derived from the averaging procedure, But- 
ler (1994) proposed instead to use long series of high quality single stations. 
This might improve the correlation between SCL and temperature. 

He showed that this was the case when using temperature series obtained 
at the Armagh Observatory in Northern Ireland 1844-1990. Since the Ar- 
magh series correlated strongly with the NH temperature, he concluded that 
this indicates that solar activity, or something closely related to it, has had a 
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dominant influence on the temperature of tlie lower atmosphere in the North- 
ern Hemisphere over the past 149 years (Butler, 1994). This investigation 
was later expanded to the period 1795-1995 by Butler & Johnston (1994), 
who found a relation tARM = 14.42 — 0.5Lsc, where Lsc is the smoothed 
length of the solar cycle determined by the 1-2-2-2-1 filter. They concluded 
that the good fit over nearly 200 years indicated that solar activity had been 
the dominant factor for nearly two centuries. 

Butler & Johnston (1996), studying the same dataset, noticed that there 
seemed to be about one solar cycle delay between the shortest cycle lengths 
and the temperature peaks. They therefore compared the smoothed value 
of SCL with the temperature 11 years later, whereby the correlations im- 
proved. They also compared the correlations with the raw, (unsmoothed) 
times of minima and maxima, filter 1-2-1, and the 1-2-2-2-1 filter, and found 
that relations were approximately the same, but the correlations improved 
significantly with the length of the filter. They interpreted this as the link 
between the solar dynamo and the mean air temperature at Armagh is a 
gradual process that becomes more evident when SCL is smoothed over sev- 
eral sunspot cycles. They also found that temperatures shifted 1 1 years back 
in time, correlated better with SCL measured between minima than between 
maxima. 

In figure 1, the length of sunspot cycles between successive minima is 
plotted inverted. Peaks in this curve preceded peaks in the Armagh temper- 
ature curve, and bumps preceded colder periods. A straight line with slope 
(3 = —(0.012 ± 0.004)yr~^ represents the average change in sunspot cycle 
length in the period 1843-1996. 

Another study by Wilson (1998) presented extremely good correlations 
between the temperature series at Armagh and the NH-temperatures (r=0.88) 
and global temperatures (r=0.83), and explained this with the influence of 
the ocean. He also found an extremely good negative correlation between the 
length of the sunspot Hale cycle of about 23 years and the Armagh temper- 
ature average. In addition to this he also found a good correlation between 
even and odd sunspot cycles and predicted t=9.24±0.47°C (90% confidence 
interval) for SC23. The observed value became 10.1°C - somewhat higher 
than predicted. 

Archibald (2008) was the first to reahze that the length of the previous 
sunspot cycle (PSCL) has a predictive power for the temperature in the next 
sunspot cycle, if the raw (unsmoothed) value for the SCL is used. Based on 
the estimated length of SC23 then being 12.6 years, considerably longer than 
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SC22 of 9.6 years^, he predicted cooling during the coming SC24 for certain 
locations. 

He demonstrated this based on a long series from de Bilt in the Nether- 
lands 1705-2000 which showed a decrease of 0.6°C per year PSCL and for 
Hanover NH with a slope of -0.73°C yr~^. Other long series he investigated 
was Portland, ME (slope: -0.70°C yr"^), Providence, RI (slope: -0.62 °C 
yr~^), and Archangel, RU with a slope -0.6 °C yr~^ PSCL. In his analysis 
he used both the times between maxima and between minima in the same 
relations (Archibald, 2010). 

In the following we will compare the relations between raw (unsmoothed) 
SCL values and temperatures in the same or in the next sunspot cycles. 
The use of unsmoothed SLC values may give a better response to the large 
change in SCL that took place from SC22 to SC23. We will also investigate 
the relations using different time lags. Our main goal is to investigate the 
use of PSCL as a tool for temperature predictions at certain locations. In 
addition we will use the relation between PSCL and temperature to estimate 
how much of the temperature variations in the series investigated, may be 
attributed to solar activity. We will also use the PSCL-relations to predict 
the future temperature changes from SC23 to SC24 - a prediction that may be 
falsified in less than 10 years. In section 2 we present data sets and methods. 
In section 3 results, and then discussion and conclusions in sections 4 and 5. 

2. Data sets and method 

Mean annual air temperatures during sunspot cycles, as far as possible 
back to 1856, which was the start of SCIO, are calculated for a selection of 
Norwegian weather stations, which have long records. Both costal and inland 
stations are represented. To avoid urban heat effects, only temperature series 
from places with small populations or stations at lighthouses are analyzed. 
These are then correlated with the length of the same (SCL) and the previous 
(PSCL) sunspot cycles. Based on relations determined between PSCL and 
temperature, temperature forecast for solar cycle 24 are calculated with 95% 
confidence interval when possible. To investigate the validity of the relations 
found, also North Atlantic temperature series from Armagh, Archangel, The 



^The recent version of SCL numbers obtained from NGDC (2011), used in this work, 
shows SCL22 = 10.0 and SCL23 = 12.2 years between minima. 
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Figure 2: The location of meteorological stations with temperature series analyzed. 

Faroe Islands, Iceland, Svalbard and Greenland are analyzed. The location 
of the stations used is shown in figure 2. 

We have also determined SCL- and PSCL-relations for composite tem- 
perature averages for Norway and we have constructed a European mean 
temperature anomaly based on 60 stations (EuropeGO), listed and shown on 
a map (figure A. 20) in the appendix. This composite series has an increasing 
number of stations included with time, as shown in figure A. 21. These sta- 
tions were selected to 1) obtain long data series, 2) obtain good geographical 
coverage, and 3) to minimize potential effects from urban heat islands. In ad- 
dition we have investigated the same relations for HadCRUTSN temperature 
anomalies. 

As an estimate of a possible solar effect related to PSCL, the coefficient 
of determination is calculated, where ris the correlation coefficient. In 
addition we have investigated present autocorrelations in the residuals when 
the relations between PSCL and temperature is employed, to examine the 
level of presence of other possible regressors. Finally we have also tested our 
forecasts by removing the last observation and made a forecast based on the 
remaining data points. 
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Figure 3: Correlations between SCL and 11 years average temperatures as function of 
time lag after the mid time of the sunspot cycle. 



3. Results 

3.1. Correlations with time lag 

As a background for the investigation of possible relations between SCL 
and temperature in sunspot periods, we determined the correlation between 
SCL and temperature for variable lags of an 11-year time window. We cal- 
culated 11 year running mean temperatures for the selected data sets, and 
correlated this with the SCL with lags from zero to 13 years, calculated from 
the middle time for each solar cycle. The starting point could also have been 
chosen as the year of solar maximum or the end year of the sunspot cycle. 
However, we selected the middle time, since this gave the possibility to check 
correlations with the same sunspot cycle (lag=0). 

A lag is defined as later in time (in the future). As input temperature 
files we used the temperature series shown in figures 4-19, lower panels. 

The result is displayed in figure 3. The North Atlantic locations (panel 
A) show a systematic better negative correlation up to a maximum negative 
for lags 9-12 years. For the 3 locations in southern Norway (panel B), the 
tendency is the same, with a maximum negative correlation between 9 and 
11 years. For Northern Norway and Archangel, the situation is somewhat 
different: The relations are not as steep as for the North Atlantic stations. 
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Three of the stations (Bod0, Troms0, and Karasjok) shows minima with 9- 
10 years lags, but a station further east, (Vard0, 31E) shows almost equal 
minima after 3 and ^10 years lags. Finally Archangel (41E), shows only one 
minimum after 2-3 years. We interpret this as the North Atlantic climate 
is dominating in Norway, except for the easternmost station (Vard0) where 
another climate zone (boreal) appears. This type of climate dominates for 
Archangel further east. This means that if we calculate temperature averages 
for stations from one geographical area, we may mix climate zones, and for 
hemispherical or global temperatures, we may have a mixture of different 
time lags and the maximum correlation may therefore not be well defined. 

Correlations with time lags for area- averaged temperatures are displayed 
in panel D, where results for Norway, EuropeGO and HadCRUT3N series are 
shown. For Norway, the maximum negative correlation is for 11 years time 
lag and well defined, while the EuropeGO average shows two minima of ^2 
and 10 years time lags, signaling the inclusion of eastern (boreal) climate 
stations. For this data set the highest correlation is after 2 years, which is 
mostly inside the same solar cycle. The HadCRUT3N series shows only one 
shallow minimum at G - 12 years time lag - apparently a mix of many time 
lags. 

3.2. Temperature trends 

The temperature series are analyzed with least square fits to the linear 
relation y = I3x -\- a, where y is temperature or temperature anomaly (for 
EuropeGO and HadCRUT3N), and x is time or length of sunspot cycles (SCL 
or PSCL). Results for data sets analyzed are shown in table 1. 

Pi is the secular temperature trend for temperatures calculated from the 
average temperature in sunspot cycles including cycle 23 which terminated 
in 2008. The starting year of the temperature series are given in column 3. 

PscL is the trend in the relation between the length of a solar cycle and 
the temperature in the same cycle. For all temperature series investigated, 
no trend significantly different from zero on the 2a level has been found for 

PsCL- 

fSpscL is the trend between the length of one solar cycle and the tem- 
perature in the next cycle. We name this the previous solar cycle length 
(PSCL) relation, since we investigate the use of the length of a solar cycle as 
predictor for the temperature in the next. 

^P5CL estimated coefficient of determination based on the PSCL 

model. 
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Table 1. Temperature trends and predictions for solar cycle 24 



Place 


Location 


Start 


Pi 






■2 

PSCL 




*c 


Conf int. 






yr 


lOOyr-l 








°C 




95% 


Oks0y 


58N,8E 


1870 


0.44±0.17 


0.05±0.14 


-0.28±0.15 


0.25 


8.4 








Utsira 


59.5N,5E 


1868 


0.75±0.21 


0.03±0.13 


-0.37±0.12 


0.45 


8.3 


6.9 


6.1:7.7 


Dombas 


62.4N,10E 


1865 


0.95±0.23 


-0.01±0.16 


-0.45±0.15 


0.55 


2.6 


0.9 


-0.1 


:1.8 


Bod0 


67.5N,ir>E 


1868 


0.78±0.20 


-0.07±0.15 


-0.42±0.11 


0.56 




3.9 


3.2 


4.7 


Troms0 


69.7,19E 


1868 


0.36±0.18 


-0.06±0.10 


-0.29±0.08 


0.56 


3.3 


2.3 


1.8 


2.8 


Karasjok 


69.5N,25.5E 


1876 


0.59±0.43 


-0.11±0.21 


-0.45±0.18 


0.43 


-1.1 


-2.6 


-3.8 


-1.4 


Vard0 


70.4N,31E 


1856 


0.91±0.18 


-0.16±0.17 


-0.38±0.12 


0.44 


2.3 


0.8 


-0.1 


:1.7 


Norway 




1900 


0.56±0.36 


0.11±0.15 


-0.30±0.12 


0.42 


1.82 


0.72 


-0.05:1.49 


Armagh 


54.3N,6W 


1867 


0.44±0.25 


0.09±0.13 


-0.29±0.13 


0.33 


10.1 


8.8 


8.0 


9.7 


Archangel 


64.BN,41B 


1881 


1.4 ± 0.3 


-0.18±0.24 


-0.51±0.21 


0.38 


1.6 


-0.6 


-1.9:0.9 


Torshavn 


62N,7W 


1890 


0.58±0.24 


-0.02±0.12 


-0.34±0.08 


0.67 


7.0 


6.9 


5.5 


6.5 


Akureyri 


65.4,18W 


1882 


1.16±0.44 


-0.06±0.25 


-0.71±0.14 


0.72 


4.2 


2.3 


1.3 


3.3 


Svalbard 


78.2,15.5E 


1914 


1.64 ±1.1 


0.69±0.41 


-1.09±0.31 


0.63 


-4.3 


-7.8 


-9.6 


-6.0 


Nuuk 


64.1, 51W 


1881 


0.7±0.6 


-0.12±0.30 


-0.65±0.24 


0.41 


-0.7 


-2.3 


-4.0 


-0.6 


Europc6() 




1856 


0.62±0.18 


-0.14±0.12 


-0.29±0.10 


0.39 


0.97 


-0.20 


-0.49:0.10 


HadCRIJT3N 




1856 


0.47±0.09 


-0.02±0.09 


-0.21±0.06 


0.49 


0.49 


-0.38 


-0.84 


0.09* 



Estimated standard deviations are given. * Uncertain because of significant correlations in residuals (see table 2) 



Usually the degrees of freedom applied in the regression analysis are equal 
to the number of observations minus the number of parameters used in the 
modeling. However, since we have investigated the correlation as function 
of time lag after the middle time of a solar cycle, and found that maximum 
correlation exist about one sunspot cycle later, the number of degrees of 
freedom is reduced by one. 

Our suggested model is rather simple: The air temperature in a sunspot 
cycle is a linear function of the length of the previous sunspot cycle (PSCL) 
in short: the t{PSCL) -model. 

Table 1 shows that for all 16 series examined, except for Oks0y, Ppscl 
is significant negative on the 95% level on the condition that the regression 
model applied gives independent residuals. 

Analysis has been performed on the residuals, which are the differences 
between the temperature observations and the fitted temperatures. To check 
possible autocorrelation in the residuals a Durbin- Watson (DW) test (Durbin 
& Watson, 1950, 1951a,b) has been applied. The results of the DW test are 
shown in Table 2. On the 5% significance level this test showed no significant 
autocorrelations in the residuals for 10 of the 16 series. Exceptions were 
the HadCRUTSN series, where significant autocorrelations were identified, 
and the Oks0y, Dombas, Bod0, Norway and Archangel series where either 
significant autocorrelations or no significant autocorrelations were identified 
(indifferent cases). 

For each of the 10 series with no autocorrelation in the residuals, the 
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Table 2. Durbin-Watson test on autocorrelations in the model residuals 

Degrees 0.05 Significance level DW test No correlation 



Series DW freedom D(L) D(U) result in residuals 



Oks0y 


1.27 


10 


0.971 


1.331 


D(L)<DW < D(U) 


Indifferent 


Utsira 


1.38 


10 


0.971 


1.331 


DW > D(U) 


Significant 


Dombas 


1.03 


10 


0.971 


1.331 


D(L) < DW < D(U) 


Indifferent 


Bod0 


1.14 


10 


0.971 


1.331 


D(L) < DW < D(U) 


Indifferent 


Troms0 


1.49 


10 


0.971 


1.331 


DW > D(U) 


Significant 


Karasjok 


1.54 


9 


0.927 


1.324 


DW > D(U) 


Significant 


Vard0 


1.37 


11 


1.010 


1.34 


DW > D(U) 


Significant 


Norway 


1.205 


7 


0.824 


1.320 


D(L)<DW<D(U) 


Indifferent 


Armagh 


1.55 


10 


0.971 


1.331 


DW > D{U) 


Significant 


Archangel 


1.24 


10 


0.927 


1.324 


D(L)<DW<D{U) 


Indifferent 


Torshavn 


1.74 


8 


0.879 


1.320 


DW > D{U) 


Significant 


Akureyri 


1.59 


9 


0.927 


1.324 


DW > D(U) 


Signiheant 


Svalbard 


2.15 


6 


0.763 


1.332 


DW < 4- D{U) 


Significant 


Nuuk 


1.67 


9 


0.927 


1.324 


DW > D{U) 


Significant 


Europe60 


1.48 


11 


1.01 


1.340 


DW > D(U) 


Significant 


HadCRUT3N 


0.86 


11 


1.01 


1.340 


DW < D{L) 


Not significant 


The Durbin-Watson lower D(L) and upper 


bounds 


D(U) numbers 


are given in Savin Sz Whit 


e (1977). 


The number of parameters k, 


is 1 in all cas 


ies since 


the regression constant is excluded. 





The number of observations in the DW test, is in all cases equal the number of degrees of freedom plus 2. 



estimated trends are significant different from since the absolute values 
of the estimated trends are greater than the estimated standard deviation, 
a, multiplied by to.025(/) - the 0.025 percentile in the Student distribution 
with / degrees of freedom. Since the residuals arc not autocorrclatcd, the 
degrees of freedom for the Student distribution are equal to the number of 
observations minus 3. As mentioned above one additional degree of freedom 
is subtracted because of inspection of the data on beforehand. Hence the 
degrees of freedom vary between 6 and 11, and the 0.025 percentile in the 
Student distribution varies between 2.20 and 2.45. 

The number of average temperature observations for Bod0 and Dombas is 
13. Then to.o25(10) = 2.23. Table 1 shows that the absolute value of the trend 
estimate for Bod0 is ~ 3.8a and for Dombas ~ 3a. The 0.025 percentile for 
the Student distribution with 3 degrees of freedom, to.025(3) = 3.18. It is 
no reason to believe that possible, but not significant autocorrelations in the 
residuals for Bod0 or Dombas will reduce the degrees of freedom from 10 to 
3 or less. Hence, the temperature trends for Bod0 and Dombas are assumed 
to be significantly different from 0. For Archangel we find Ppscl/c^iSpscl ~ 
2.43 > to. 025(10) = 2.23. Hence, also the Archangel temperature model is 
assumed to give significant results on a 95 % significance level. 

Analysis of the Oks0y series does not give satisfactory results. Hence, 
because of a possible dependency in the residuals and a significance level 
lower than for the other series, it is not possible to estimate a 95 % confidence 
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TromsB: Temperature corrected for sunspot cycle length effect 




Figure 4: Troms0, Norway (Coastal city 70N): Lower panel: Annual mean temperatures 
and mean temperatures in sunspot cycles (dots) with linear trend. The broken line shows 
temperatures measured so far in SC24. The center two panels show temperatures in 
sunspot cycles versus SCL (left) or PSCL (right). In the right center panel the 95% 
confidence interval is indicated with dotted lines. The lower panel shows the forecasted 
value for SC24 (diamond) with standard deviation error bar. The top panel shows the 
residuals from the t{PSCL) model, and represents what may be left for other regressors 
and noise to explain. 

level for the forecasts in this case. 

For the average temperature of Norway we get PpsclIctppscl — 2-50 > 
^0.025(7) = 2.377, which indicates that this is a significant result on the 95% 
significance level. 

The absolute value of the estimated HadCRUT3N temperature trend is 
~ 3.5(7. The number of temperature observations is 14 and the number of 
degrees of freedom 11. Then ^0.025(11) = 2.20. Even with a reduction of 
number of freedoms from 11 to 3, the temperature trend will be significant 
different from 0. However, in this case there is reason to be more careful 
since significant autocorrelations in the residuals are identified. 

The regression analysis indicates clearly that the temperatures in the 
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DombSs: Temperature corrected for sunspot cycle length effect 
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Figure 5: Dombas, Norway (Inland, mountain valley 62N), explanation as in figure 4 



PSCL model have significant trends different from for 14 out of the 16 data 
series. It is also reasons to believe that the trend is significantly different from 
for the HadCRUT3N temperature, but there is some uncertainty because 
of significant autocorrelations in the residuals. 

A 95% confidence interval of the estimated trend is given by respectively 
subtracting and adding the factor a x to.025(/) to the estimated trend. This 
is the 95% confidence interval for the temperature forecasts shown in the 
last column in table 1. For Oks0y no confidence interval is estimated. The 
next 14 confidence intervals are considered to be precise, while the last one 
for HadCRUT3N can be considered to be an approximate 95% confidence 
interval. We conclude that 15 of our 16 series support the PSCL regression 
relation. 

The coefficient of determination (r^) indicates how much of the variance 
can be attributed to the regressor. Our results show that it varies from 0.25 
at Oks0y, at the southern tip of Norway, to 0.72 at Akureyri at the bottom 
of a large fjord in Northern Iceland. The highest values are found on islands 
and in northwestern coastal regions of the North Atlantic. Lower values are 
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Norway: Temperature corrected for sunspot cycle length effect 
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Figure 6: Norway, average yearly temperatures, explanation as in figure 4. 



found inland, and distant from the main Atlantic Ocean currents. 

3.3. Local temperature series 

The results for two stations in Norway, one coastal north (Troms0) and 
one inland south (Dombas), are shown in figures 4 and 5. For these two 
locations, and the others investigated and shown as figures 8-18 the Pscl- 
values are not significantly different from zero. This is demonstrated in 
the center left panels in the figures, which show the relations between the 
temperature and length of the same sunspot cycle. On the other hand, 
significant relations are found between PSCL and temperatures, which are 
shown in the center right panels. The residuals from the PSCL relations 
(sorted time wise) shown in the top panels, indicate the variability that 
remains to be explained either by other regressors or noise. 

For Troms0 (figure 4) the secular temperature trend before apphcation 
of the PSCL-model is reduced from /3i = 0.0036°Cyr-i to P2 = 0.00096 ± 
0.0004°Cyr^^, where is the secular trend in the residuals after removal 
of the sunspot cycle effect (PSCL-model). The coefficient of determination 



15 



EuropeSO: Temperature corrected for sunspot cycle length effect 
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Figure 7: Temperature anomaly for 60 stations on land in Europe, mostly located outside 
large cities. A list of the stations is given in the Appendix, explanation as in figure 4. 

= 0.56 indicates that more than half of the temperature variations may 
be attributed to solar activity. 

For Dombas, situated in a mountain valley in southern Norway, (figure 
5) the coefficient of determination is = 0.46, indicating that slightly less 
than half the temperature variations are related to solar activity as modeled 
by PSCL. For Utsira, which is a hghthouse on the western coast of Norway 
(figure 9) = 0.45. The similar results indicate that the Atlantic currents 
may infiuence both locations. On the other hand, for Oks0y (see figure 8), a 
lighthouse on the southern tip of Norway, we found that the PSCL-relation 
was not significant on the 95% level, and this is refiected in = 0.25, which 
is the smallest solar contribution for the locations investigated. This location 
is influenced by different ocean streams that may cancel out the effect. 

In figures 4-19 are also shown (with diamonds) predicted average tem- 
peratures for SC24 based on the PSCL-model. Temperatures measured for 
2009 and 2010 are shown as broken lines. For many locations the mean an- 
nual temperatures have decreased in the first years of SC24. In some of the 
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Oks0y: Temperature corrected for sunspot cycle length effect 
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Figure 8: Oks0y lighthouse, Norway, average yearly temperatures, explanation as in figure 
4. The confidence interval indicated in the center right panel is probably less than 95% in 
this case. 



figures yearly observations done before the first SC included in the analyses 
are connected with broken lines. 

The average solar (PSCL) contribution to the variability for the 4 north- 
ern Norwegian stations in table 1 is 50%, compared to 42% for the 3 stations 
located in southern Norway. Moving further west and north, we find that 
the stations at Torshavn, Akureyri and Svalbard (figures 15-17), have the 
maximum solar contribution with = 0.63 — 0.72. 

For Archangel (figure 13) we get Ppscl = —0.51 ± 0.21, which is in line 
with the result of Archibald (2010) who got (5 = —0.6, based on PSCL calcu- 
lated from both the time between maxima and between minima. Archangel 
displayed the highest correlation with a time lag of only 2 years (figure 3). 
This may indicate a boreal climate different from the Atlantic, because of 
the shorter reaction time to a solar pulse. 

For Armagh (figure 14) we get (5pscL — —0.29 ± 0.13, which is smaller 
than the value /3 = —0.5 determined by Butler & Johnston (1994). Their 
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Utsira: Temperature corrected for sunspot cycle length effect 
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Figure 9: Utsira lighthouse, Norway, average yearly temperatures, explanation as in figure 
4. 
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Bodo: Temperature corrected for sunspot cycle length effect 
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Figure 10: Bod0, Norway, average yearly temperatures, explanation as in figure 4. 
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Karasjok: Temperature corrected for sunspot cycle length effect 
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Figure 11: Karasjok, Norway, average yearly temperatures, explanation as in figure 4. 



larger value may come from their smaller range in SCL values due to the 
1-2-2-2-1 smoothing applied by them. 

The secular trend is determined with least significance for Svalbard (78N), 
which also has the largest trend, due to Arctic amphfication. The Svalbard 
temperature record is based on measurements at the coast of the large fjord 
Isfjorden, and the temperature is therefore affected by both regional and local 
sea ice conditions. It is somewhat surprising that the coefficient of determi- 
nation = 0.63, which indicates that solar activity plays an important role 
for this Arctic location far North of the Polar Circle. The series for Svalbard 
is also much shorter than the other series, but is still the longest in the high 
Arctic. The temperature measured at Svalbard has already shown sign of 
decline as predicted. 

A more detailed analysis of the Svalbard temperature series has been per- 
formed by Solheim et al. (2012), who repeated this analysis with weighted 
averages and also analyzed each season separately. They confirmed the val- 
ues of the correlation coefficient r with a bootstrap analysis, and got similar 
results as reported here for the coefficient of determination r^. Only the 
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Varde: Temperature corrected for sunspot cycle length effect 
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Figure 12: Vard0, Norway, average yearly temperatures, explanation as in figure 4. 

yearly average and the winter season (DJF) temperatures showed no auto- 
correlations in the residuals. The winter temperatures gave = 0.67, which 
indicates that the solar influence takes place also when the Sun is completely 
under the horizon. Indirect lagged contributions by means of the ocean cur- 
rent and warm air advected from the south are the only possible explanations 
for the correlations determined. 

For Nuuk on Greenland (figure 18), the secular trend is not signifi- 
cant, and disappears almost completely (^2 — 0.0006 ± 0.0002°Cyr~^) when 
corrected for solar activity. The Nuuk station is located near the coast in 
West Greenland, and is affected by regional sea surface temperatures. The 
sea between the Labrador island and West Greenland has shown a marked 
warm anomaly for most of 2010. This has resulted in a temperature increase 
the last two years (figure 18), opposite the prediction in table 1. 

3.4- Area averaged air temperature series 

The yearly average temperature for Norway calculated from 1900 on- 
wards, is shown in the bottom panel of figure 6. A significant relation with 
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Archangel: Temperature corrected for sunspot cycle length effect 
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Figure 13: Archangel, Russia, average yearly temperatures, explanation as in figure 4. 
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Armagh: Temperature corrected for sunspot cycle length effect 
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gure 14: Armagh, Northern Ireland, average yearly temperatures, explanation as 
ure 4. 
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Torshavn: Temperature corrected for sunspot cycle length effect 
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Figure 15: Torshavn, Faroe Islands, average yearly temperatures, explanation as in figure 
4. 
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Akureyri: Temperature corrected for sunspot cycle length effect 
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Figure 16: Akureyri, Iceland, average yearly temperatures, explanation as in figure 4. 
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Svalbard: Temperature corrected for sunspot cycle length effect 
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Figure 17: Svalbard (Longyearbyen), average yearly temperatures, explanation as in figure 
4. 
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Nuuk: Temperature corrected for sunspot cycle lengtti effect 
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Figure 18: Nuuk, Greenland, average yearly temperatures, explanation as in figure 4. 
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HadCRUT3N: Temperature corrected for sunspot cycle length effect 
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Figure 19: HadCRUTSN, average yearly temperature anomalies, explanation as in figure 
4. 



PpscL — —0.30 ± 0.12 reduces the secular trend with a large fraction to 
/32 = 0.0020 ± 0.0002°Cyr-^ Solar activity may contribute at least 42%. 

For the average temperature for 60 European stations, mostly located 
outside big cities, we find a significant relation with (ipscL = —0.29 ± 0.10 
(figure 7) and = 0.39. HadCRUT3N temperature anomalies gives Ppscl = 
—0.21 ± 0.06, and = 0.49, but in this case also other regressors may 
contribute significantly. 

As seen in figures 6 and 7, the Norwegian and Europc60 average tem- 
peratures have already started to decline towards the predicted SC24 values, 
while the HadCRUT3N temperature anomaly has shown no such decline yet 
(figure 19). 

3. 5. Testing the PSCL-model forecast 

To check the robustness of the forecast with the PSCL-model we have 
for all 16 datasets removed the last entry, which is the SC23 temperature, 
and used the remaining observations for generating a forecast of t(SC23). 
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The result is that the forecast for 11 of the 15 series with 95% confidence 
intervals are within this interval. For 3 of the series (Dombas, Svalbard 
and HadCRUTSN, the forecasted temperature is 0.1°C above the forecasted 
range, and in the case of Bod0, the deviation is 0.2°C. If the SC23 had been 
a fraction shorter, all observed values would be inside the 95% range. It 
should be remarked that SC22 for a long time was hsted in NGDC (2011) 
with a length of 9.6 or 9.7 yr, which was used in other analysis (Archibald, 
2008; Thejll, 2009). According to earher lists (Thejll, 2009) May 1996, which 
is the mathematical minimum for SC93, was in the beginning used, while 
a consensus among solar scientists moved the minimum to September 1996, 
which made SC22 longer (10.0 yr) and SC23 shorter (12.2 yr). This difference 
in definition of SCL from the original by Waldmeier (1939) may explain the 
observations outside the confidence range. 

4. Discussion 

The correlations found between the average temperature in a solar cycle 
and the length of the previous cycle, indicates a possible relation between 
solar activity and surface air temperature for the locations and areas inves- 
tigated. 

The weak or non significant, correlations between the temperature and 
length of the same cycle, and the much stronger correlations for time lags of 
the order 10-12 years, makes the length of a solar cycle a good predictor for 
the average temperature in the next cycle. 

The smoothing filters 1-2-2-2-1 and 1-2-1 or a length jump (sudden in- 
crease of the SCL), may be the reason that the previous good correlations 
between solar cycle length and Northern hemisphere land-temperature dis- 
appeared after 2000 (Thejll & Lassen, 2000; Thejll, 2009). 

We expect that the optimum delay is not exactly one solar cycle at some 
locations on the Earth, presumably due to oceanographic effects. This is 
confirmed by our determination of the correlations between SCL and tem- 
perature lags which appears to be at a maximum negative value for 10-12 
year lags around the North Atlantic, but shorter at eastern locations as Vard0 
on the Barents sea and Archangel in Northern Russia (figure 3). 

Variations in temperature lags may not give precisely one solar cycle lag 
as the optimum relations between area averaged temperatures and SCL, but 
both for the Europe60 station average and HadCRUT N acceptable correla- 
tions are found with PSCL - and reasonable good predictions are possible. 
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The PSCL- model predicts a temperature drop of 1.3 to 1.7°C for single 
Norwegian stations analyzed from SC23 to SC24. For the average Norwegian 
and Europa60 temperatures the temperature drop from SC23 to SC24 is 
1.1-1.2 °C. For HadCRUT3N the predicted temperature drop is 0.9°C. 95% 
confidence intervals for the predicted temperatures in SC24 are given in table 
1, last column. The Arctic cooling as predicted here may be converted into 
a global cooling, which is a factor 2-3 lower due to the Artie amplification of 
temperature differences (Moritz et al., 2002). Tnis means a global cooling of 
the order 0.3-0.5°C. We may also expect a more direct cooling near Equator 
due to the response to reduced TSI with the weaker solar cycles in the near 
future (Perry, 2007; Mehl et al, 2009; Richards et al., 2009). 

The PSCL relation is determined for the period 1850-2008 when the PSCL 
on average has shortened from cycle to cycle in relative small steps, and the 
Earth has warmed. The large increase in SCL from SC22 to SC23 signals a 
temperature drop, which may not come as fast as predicted because of the 
thermal inertia of the oceans. The warming has taken place over 150 years - 
cooling of the same order may require some decades to be realized. 

The temperatures in the North Atlantic and on adjacent land 
controlled by the Atlantic Multidecadal Oscillation (AMO), which has in 
the instrumental period from 1856 exhibited a 65-80 yr cycle (0.4°C range), 
with warm phases at roughly 1860-1880 and 1930-1960 and cool phases dur- 
ing 1905-1925 and 1970-1990 (Kerr, 2000; Gray et al., 2004). This period 
is proposed related to a 74 yr period, which is a sub harmonic of the 18.6 
yrs lunar-nodal-tide (Yndestad et al., 2008). The AMO may affect the tem- 
peratures at the locations we have analyzed, but a preliminary investigation 
by the authors have found considerable variations in dominant periods from 
place to place, which we interpret as due to local variations related to ge- 
ography, dominant weather patterns, and distance from the ocean currents. 
Inclusion of periodic variations may improve the forecasts, and we plan to 
investigate this further. 

The coefficient of determination is a measure of the relative contribu- 
tion of solar activity modeled by the PSCL-model, to the temperature varia- 
tions. The range of is from 0.25 (Oks0y) to 0.72 (Akureyri). The highest 
values (Torshavn: 0.62, Akureyri: 0.72 and Svalbard: 0.63) are found in or 
close to the North Atlantic ocean currents. This points to ocean currents as 
the mechanism of transport of the heat generated at southern locations by 
solar radiation. 

Even if the correlation is rather high between the actual temperature and 
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the fitted temperature by use of regression modeling, tlie question arises of 
possibilities to improve the models. Other regressors, which may be included, 
are atmospheric and oceanic dynamics unrelated to the Sun or with different 
time lags, volcanic activity, GHG forcing, and anthropogenic activity of var- 
ious kinds. The contribution from other regressors should be determined for 
each location, and may be quite different from area-averaged regressors. This 
is a question left for future investigations. However, there are clear limita- 
tions for including additional explanatory variables simply because there are 
limited numbers of observations available. Additional variables will improve 
the model fitting, but reduce the degree of freedom. 

This analysis shows significant dependency between the previous sunspot 
cycle length and the temperature. The established model is able to make 
significant forecasts with 95% confidence limits for the present sunspot cycle. 
There are reasons to believe that these results could be fundamental in further 
development of long-term forecasting models for the temperature. 

Looking at figure 1, which shows the variation in the length of solar cycles, 
we realize that short cycles like the one that ended in 1996, have only been 
observed 3 times in 300 years. After the shortest cycles, sudden changes 
too much longer cycles have always taken place, and thereafter there is a 
slow shortening of the next cycles, which take many cycles to reach a new 
minimum. This recurrent pattern tells us that we can expect several long 
cycles in the next decades. Analysis of the SCL back to 1600 has shown a 
periodic behavior with period 188 yr, now entering a phase with increasing 
SCL the next ~75 yr (Richards et al., 2009). 

de Jager & Duhau (2011) concludes that the solar activity is presently 
going through a brief transition period (2000-2014), which will be followed 
by a Grand Minimum of the Maunder type, most probably starting in the 
twenties of the present century. Another prediction, based on reduced solar 
irradiance due to reduced solar radius, is a series of lower solar activity cycles 
leading to a Maunder like minimum starting around 2040 (Abdussamatov, 
2007). 

A physical explanation for the correlations between solar activity parame- 
ters as the SCL and the temperatures on the Earth, and a possible decoupling 
of this relation the last few years has been investigated by Solanki & Krivova 
(2003). Their conclusion is that even if the proxies and direct measure- 
ments are scaled such that statistically the solar contribution to climate is as 
large as possible in the period 1856-1970, the Sun cannot have contributed 
to more than 30% of the global temperature increase taken place since then. 
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irrespective of which of the three considered channels is the dominant one de- 
termining Sun-cUmate interactions: Tropospheric heating caused by changes 
in total solar irradiance, stratospheric chemistry influenced by changes in the 
solar UV spectrum, or cloud coverage affected by the cosmic ray flux. 

Prom correlation studies of 7 (not all global) temperature series for the 
period 1610-1970 de Jager et al. (2010) found a solar contribution of 41% 
to the secular temperature increase. Our results are somewhat higher for 
Northern Hemisphere locations in the period 1850-2008. 

Analyzing global temperature curves for periodic oscillations Scafetta 
(2010) concludes that the chmate is forced by astronomical oscillations re- 
lated to the Sun, and at least 60% of the warming since 1970 can be related 
to astronomical oscillations. Looking at our figures 4-19, we can see sign of 
quasi-periodic variations in the lower panels, which have more or less disap- 
peared in the upper panels. We may therefore suggest that SCL in some way 
is related to astronomical forcing. 

Satellite observations by the Spectral Imager Monitor (SIM) indicate 
that variations in solar ultraviolet radiation may be larger than previously 
thought, and in particular, much lower during the recent long solar minimum. 
Based on these observations Ineson et al. (2011) have driven an ocean-climate 
model with UV irradiance. They demonstrate the existence of a solar chmate 
signal that affects the NAO (North Atlantic Oscillation) and produced the 
three last cold winters in Northern Europe and in the United States. 

5. Conclusions 

Significant linear relations are found between the average air temperature 
in a solar cycle and the length of the previous solar cycle (PSCL) for 12 out 
of 13 meteorological stations in Norway and in the North Atlantic. For 9 
of these stations no autocorrelation on the 5% significance level was found 
in the residuals. For 4 stations the autocorrelation test was undetermined, 
but the significance of the PSCL relations allowed for 95% confidence level 
in forecasting for 3 of these stations. Significant relations are also found 
for temperatures averaged for Norway, 60 European stations temperature 
anomaly, and for the HadCRUT3N temperature anomaly. Temperatures 
for Norway and the average of 60 European stations showed indifferent or no 
autocorrelations in the residuals. The HadCRUTSN series showed significant 
autocorrelations in the residuals. 
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For the average temperatures of Norway and the 60 European stations, 
the solar contribution to the temperature variations in the period investigated 
is of the order 40%. An even higher contribution (63-72%) is found for 
stations at Faroe Islands, Iceland and Svalbard. This is higher than the 7% 
attributed to the Sun for the global temperature rise in AR4 (IPCC, 2007). 
About 50% of the HadCRUTSN temperature variations since 1850 may be 
attributed solar activity. However, this conclusion is more uncertain because 
of the strong autocorrelations found in the residuals. 

The significant linear relations indicate a connection between solar activ- 
ity and temperature variations for the locations and areas investigated. A 
regression forecast model based on the relation between PSCL and the aver- 
age air tempereaure is used to forecast the temperature in the newly started 
solar cycle 24. This forecast model benefits, as opposed to the majority of 
other regression models with explanatory variables, to use an explanatory 
variable - the solar cycle length - nearly without uncrtainty. Usually the ex- 
planatory variables have to be forecasted, which of cause induce significant 
additional forecasting uncertainties. 

Our forecast indicates an annual average temperature drop of 0.9°C in 
the Northern Hemisphere during solar cycle 24. For the measuring stations 
south of 75N, the temperature dechne is of the order 1.0-1.8°C and may 
already have already started. For Svalbard a temperature decline of 3.5°C 
is forecasted in solar cycle 24 for the yearly average temperature. An even 
higher temperature drop is forecasted in the winter months (Solheim et al., 
2012). 

Artie amplification due to feedbacks because of changes in snow and ice 
cover has increased the temperature north of 70N a factor 3 more than below 
60N (Moritz et al., 2002). An Artie cooling may relate to a global cooling 
in the same way, resulting in a smaller global cooling, about 0.3-0.5 °C in 
SC24. 

Our study has concentrated on an effect with lag once solar cycle in order 
to make a model for prediction. Since solar forcing on climate is present on 
many timescales, we do not claim that our result gives a complete picture of 
the Sun's forcing on our planet's climate. 
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Appendix A. The Europe60 stations 



Table 3. Location of Meteorological stations included 
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Figure A. 20: The location of the Europe60 stations 



39 



EuropeSO: No of stations 




1840 1860 1880 1900 1920 1940 1960 1980 2000 2020 

Year 



Figure A. 21: The variation in numbers of stations in the Europe60 
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